The non-Markovian stochastic Schròdinger equation for open 

systems 



Lajos Diósi 1 ' 2 '* and Walter T. Strunz 1,3,t 
1 Department of Physics, Queen Mary and Westfield College, University of London, 

Mile End Road, London El 4NS, UK 
2 KFKI Research Institute for Particle and Nuclear Physics, H-1525 Budapest 114, PO Box 49, 

Hungary 

3 Fachbereich Physik, Universitat GH Essen, D-45117 Essen, Germany 

(June 18, 1997) 

Abstract 

We present the non-Markovian generalization of the widely used stochastic 
Schròdinger equation. Our result allows to describe open quantum systems 
in terms of stochastic state vectors rather than density operators, without 
Markov approximation. Moreover, it unifies two recent independent attempts 
towards a stochastic description of non-Markovian open systems, based on 
path integrals on the one hand and coherent states on the other. The latter 
approach utilizes the analytical properties of coherent states and enables a 
microscopic interpretation of the stochastic states. The alternative first ap- 
proach is based on the general description of open systems using path integrals 

as originated by Feynman and Vernon. 
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In the last few years the description of open quantum systems in terms of stochastic 
Schròdinger equations has received remarkable attention. They are now widely used in 
different fields (measurement theory, quantum òptics, quantum chaos, sòlid states [0-0]), 
wherever quantum irreversibility matters. They do not only serve as a fruitful theoretical 
concept but also as a practical method for computations in the form of quantum trajec- 
tòries. Up to now, however, the Markov approximation was believed to be essential for a 
stochastic description ||. For systems where non-Markovian effects are inevitable, as for 
non-equilibrium relativistic fields, especially in quantum cosmology f| or sòlid state physics 
||Ï0|,|ri|l , an advantageous stochastic pure state description was missing. This Letter presents 
an exact non-Markovian stochastic Schròdinger equation. 

Traditionally, open systems are described by the reduced density operator p sys (t) = 
Tr em) (ptot{t)) , obtained from the total density operator by tracing over the environmental 
degrees of freedom. In Markov approximation, it is well known that the system dynamics 
can either be described by a master equation for the reduced density operator p sys (t), or 
alternatively, by a stochastic Schròdinger equation for state vectors \ipz(t)) [0-0] . In this 
latter approach, the reduced density operator is recovered as the ensemble average over these 
stochastic pure states: 



Here, \ipz(t)) indicates the solution of the stochastic Schròdinger equation with a particular 
realization of the - in this case Wiener - stochastic process Z(t), the mean Mz[- ■ •] refers to 
the ensemble average over these processes. The states \ipz{t)) may or may not be normalized 
depending on whether one utilizes the nonlinear || or linear [0] version of the stochastic 
Schròdinger equation. Both the linear and the non-linear equation lead to the correct reduced 
density operator according to (0) and they are mathematically equivalent by virtue of a 
redefinition of the stochastic processes Z(t) |ï^| . 

It is the aim of this paper to demonstrate that a stochastic decomposition just like (p]) also 
holds in the general case, without any approximation, in particular without Markov approx- 
imation. We derive the linear version of the relevant non-Markovian stochastic Schròdinger 
equation, the corresponding non-linear, norm preserving theory can be found similar to the 




(1) 



Markovian case |IÏ3| . 



Our result can be based on two recent independent approaches to a stochastic description 
of non-Markovian open systems |H|,P^| ■ One approach [|Ï4J] uses coherent states and has 



the advantage of offering an interpretation of the solutions of the stochastic Schròdinger 
equation from first principies. The other [TTS) is based on the Feynman-Vernon approach to 



open systems using path integrals |Ï6| and is vàlid for arbitrary temperatures. 



To be specific, we use a Standard model of open system quantum mechanics, a system 



coupled linearly via position coupling to an environment of harmònic oscillators |JT6 



H tot = H sys (q, p)-q^2 XíQí + E [P- + \ m i^Ql ] , ( 2 ) 

i i \ Zm i 1 ) 

where for later purposes we also introduce F = Ya XíQíi the force acting on the system as 
induced by the environment. We assume a factorized total initial density operator p íoí (0) = 
\'4 , o){ r 4 ) o\ ® Pt- The environment oscillators are assumed to be in a thermal initial state pV 
at temperature T, and, for simplicity, the system is assumed to be in a pure state |·0o)(V'o|· 
The time evolution of the total system is determined by the unitary von-Neumann equation 
Ptot = -i[Htot,Ptot\- 

Without any approximation, the reduced density operator p sys (t) of the model (g) can 
be represented as the ensemble average ([[]) of stochastic pure states \ifjzit)). They are the 
solutions of the following non-Markovian stochastic Schròdinger equation 

\Mt)) = -iH sys {q,p)\Mt)) +iqZ(t)\ij>z(t)) +i§ fd* <x(t,s) 6 ~j§^, (3) 

which is the main result of this paper. Equ. (||) is a stochastic equation, since it depends on 
a stochastic process Z(t) specified below. It is also non-Markovian due to a memory term 
involving the dependence of the current state \ipz(t)) on earlier noise Z(s), describing the 
(delayed) back reaction of the environment on the system. 

The dynamical properties and the temperature of the environment determine the memory 
kernel a(t,s) = J2% (xl /^ m i u 'i) [coth(huii/2kBT) cosuiit — s) — i sin(í — s)] |Ï6| . It can be 
regarded as the force correlation function a(t, s) = Tr (^F(t)F(s)pTj , where F(t) is the 
Heisenberg operator of the force of the model (0) of the undisturbed environment. This 
memory kernel also determines the probability distribution of the stochastic processes Z(t) 



entering the non-Markovian stochastic Schrodinger equation ([3]). They are coloured complex 
Gaussian processes with properties 

M z [Z(t)} = 0, M z [Z(t)Z(s)] = and M z [Z(t)Z*{s)\ = a*(t, s), (4) 

designed in such a way as to mimic the effect of the quantum force F(t). 

In the next two parts we will prové the central assertion of this paper: the solutions 
°f the non-Markovian stochastic Schrodinger equation (^) reproduce the exact re- 
duced density operator p sys (t) if one takes the ensemble average over the stochastic processes 
Z(t) according to ([]]). 

The first proof uses path integrals. The propagator J(t; 0) of the reduced density operator 
of the model (0) can be found in Feynman and Vernon's original paper |Ï6| |, 



■J(q,q',t;%,q'o,0) = J T>[q] J V[q']exp{iS sys [q] -iS sys [q']} F[q,q'\, (5) 

with the influence functional T[q, q'} encoding the effects of the environment on the system. 
It has been shown recently [|ÏJ| that the propagator (|5]) allows an exact stochastic decom- 
position using the coloured complex Gaussian stochastic processes Z(t) |17j with properties 

(i), 

J(q, q', t; q , q' Q , 0) = M z [G z (q, t; q , 0)G* z (q', t; q' , 0)] . (6) 
In |p.5|| , the stochastic propagators G z (t; 0) were given by their path integral representation 



Gz(q, t; q , 0) = J V[q\ exp iS sys [q]+iJ^ drq T Z(r) - dr do q T a(r,a) q a 



(7) 



where we concluded that the states 

\M*)) = Gz{t-M%) (8) 

recover the reduced density operator p sys {t) according to (|T]). 

Now we derive the Schrodinger equation corresponding to the stochastic propagator (^): 
first we take the time derivative on both sides of equation (^). The action S sys [q] in the 
exponent of the path integral expression (|?p leads to the system Hamilton operator H sys 
in the Schrodinger equation, the stochastic integral J drq T Z(r) adds the stochastic driving 



term qZ(t) in (^|). The only complication arises from the contribution of the double time 
integral in the exponent of the action in the path integral expression ([7|). Using functional 
differentiation, it is straightforward to show that it leads to the remaining memory term in 
(||). This completes the desired proof of the equivalence of the stochastic propagator ([?[) 
and the non-Markovian stochastic Schròdinger equation (^). With (|5[) we conclude that the 
solutions \ipz(t)) of our equation (||) do in fact recover the reduced density operator of the 
model (|^) by taking the ensemble average ([]]) over the processes Z(t). 

It is easy to show - see also |Ï5| - that our equation reduces to the well-known (linear) 



Markovian stochastic Schròdinger equation with complex Wiener noise in the limit of white 
noise. 

An alternative derivation of the non-Markovian stochastic Schròdinger equation (|3|) uses 
a coherent state basis for the environmental degrees of freedom. This route to a stochastic 



description of non-Markovian open systems was taken in [14]. We use the unnormalized 



Bargmann coherent states \a) [|T^] which are analytical in a and satisfy the completeness 
relation / d 2 a e~' a ' 2 \a)(a\ = 1, where d 2 a = dKea dlrna/n. Assume an expansion of the total 
state vector using a Bargmann basis \a) = |ai) ® \a2) ® ■ ■ ■ for the environmental degrees of 
freedom a = (d, a 2 , . . .), 

\^tot(t)} = I d 2 a e" |a|2 |^(í)> ® |a). (9) 

The states \ip a *(t)) of the system correspond to a particular 'configuration' \a) of the envi- 
ronment. Tracing over the environment, using the representat ion (||) for the total state, we 
find that the reduced density operator takes the form 



Psys(t) = Jd 2 ae- W \^(t))(^(t)\ = M a [\^(t)){^(t)\\. (10) 

For the last expression we regard the coherent state variables a as classical stochastic vari- 
ables with Gaussian distribution M a [. . .] — j d 2 a[. . .] exp{ — |a| 2 }. 

We now turn our attention to the time evolution of the total state using the coherent 
state representation @. First, we rewrite the total Hamilton operator (§) in terms of the 
creation and annihilation operators à\ and ài of the environment oscillators, and also change 



to a Heisenberg representation of the environmental part of the total Hamiltonian. Tlien we 



find triat the system part \ip a *(t)) of the total state obeys the Schròdinger equation [IB 



= -4 s ii4í))+€^(^«:+^^) ivv(í)}. (ii) 

We now restrict ourselves to the zero temperature case (T = 0), for which the initial 
condition |^ a »(0)) = \i/; ) holds for all configurat ions a. In this coherent state approach, the 



non-zero temperature case is non-trivial and will be treated elsewhere [|Ï3[ . 

In equ. flïOD the reduced density operator is expressed naturally as an ensemble aver- 
age over the pure states \ip a »(t)). Accordingly, the evolution equation (|ïl~l) represents the 
stochastic Schròdinger equation for these states. 

Remarkably, this construction is identical to the non-Markovian stochastic Schròdinger 
equation (|3|). To see the equivalence, we define stochastic processes 

(12) 

A simple calculation shows that these processes are - for zero temperature - realizations of 
the coloured complex Gaussian stochastic processes ([|): 

M a [Z a (t)} = 0, M a [Z a (t)Z a (s)} = and M a [Z a (t)Z* a (s)} = a* T=0 (t, s). (13) 

Moreover, using the chain rule we find 

£^==e-^A = / d a £^- c ^.(*-> * = í dsa T=0 (t,s)-^--. (14) 
^ ^2miüJi da* J ^ 2miUi SZ a (s) J dZ a (s) 

Replacing the expressions (0) and (p!4|) in the coherent state stochastic Schròdinger equation 
(|ïl"|), we recover the non-Markovian stochastic Schròdinger equation @, our bàsic result. 

According to (|9|), the solutions \ipz{t)) of the non-Markovian stochastic Schròdinger 
equation (|3]) represent the part of the total state which corresponds to a certain classical 
'configuration' a of the environment. The link between the process Z(t) and the configuration 
a is given by (|T2|). The interpretation of the classical 'configuration' a of the environment 
(and thus of the stochastic processes Z(t) too) can be given in the broader framework of 



hybrid densities |L9 . 



We have presented the non-Markovian stochastic Schròdinger equation. It allows the 
description of open quantum systems in terms of stochastic state vectors rather then density 
operators, without relying on the Markov approximation. 

Two alternative derivations are given. First, we establish the connection to the path 
integral approach to open systems as initiated by Feynman and Vernon. The non-Markovian 
stochastic Schròdinger equation reflects a stochastic decomposition of the propagator for the 
reduced density operator into stochastic propagators for state vectors. Secondly, we establish 
the connection to a coherent state description of the environment, allowing a microscopic 
interpretation of the stochastic states \ipz{t)). 

Our theory is exact and the model (or some straightforward generalization of it) appears 
in many areas of physics (electron-phonon interaction, the spin-boson model, the whole of 
quantum òptics, relativistic field theories as relevant as QED). These problems can now 
be phrased in the language of stochastic evolution equations in Hilbert space, without ap- 
proximation. The success of quantum trajectory methods in the Markovian case suggests 
that our result also represents a promising step towards an effective numerical algorithm 
for non-Markovian reduced dynamics. As a fundamental concept, stochastic pure state rep- 
resentations no longer depend on the Markov assumption. Starting from our exact result, 
approximations like perturbation expansions are possible, further simplifying the description 
of non-Markovian reduced dynamics. 
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